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ABSTRACT 

Magnetars have been suggested as the most promising site for the origin of observed soft gamma-ray 
repeaters (SGRs) and anomalous X-ray pulsars (AXPs). In this work we investigate the possibility 
that SGRs and AXPs might be observational evidence for a magnetic phase separation in magnetars. 
We study magnetic domain formation as a new mechanism for SGRs and AXPs in which magnetar- 
matter separates into two phases containing different flux densities. We identify the parameter space 
in matter density and magnetic field strength at which there is an instability for magnetic domain 
formation. We conclude that such instabilities will likely occur in the deep outer crust for the magnetic 
Baym, Pethick, and Sutherland (BPS) model and in the inner crust and core for magnetars described 
in relativistic Hartree theory. Moreover, we estimate that the energy released by the onset of this 
instability is comparable with the energy emitted by SGRs. 

Subject headings: gamma rays: stars — instabilities — stars: interiors — stars: magnetic field - 
stars: neutron — X-rays: stars 



1. INTRODUCTION 

Soft gamma-ray repeaters (SGRs) are compact objects 
undergoing episodic instabilities which produce super- 
Eddington X-ray outbursts. Up to now 6 SGRs (4 con- 
firmed and 2 candidates) have been observed 3 . They are 
believed to be a new class of 7-ray transients that are 
different from the source of ordinary gamm a-ray bursts. 
Obse rvations of the spin-down time-scale (|Kouveliotiul 
119981 ) have confirmed the fact that these SGRs are newly 
born neutron stars with a very large surface magnetic 
field (B ~ 10 15 G). Such star s have been named mag- 
netars (jDuncan fc Thompson! 119921 iKouveliotiu et al.l 

f }Q3T). About 10 AXPs are also catego rized as magnetars 
[aspi 20071: Ivan Paradiis et al.| [J995). Even though the 
magnetar model is generally accepted as the paradigm 
for both SGRs and AXPs, it is not easy to explain both 
objects simulta neousl y by a consistent set of parameters. 
Woods et al| (|2001l ) have reported evidence for a sud- 
den magnetic field reconfiguration in SGR 1900+14 dur- 
ing the giant flare of August 27, 1998. This scenario 
requires a reorganization of the magnetic field both in- 
side and outside the star. Sharp field gradients are pos- 
tulated to create a fracture of the ri gid outer crust of 
the neutron star. Cheng et al| (|1995| ) have shown that 
SGR events and earthquakes share four distinctive sta- 
tistical properties: 1) power-law energy distributions; 
2) log-symmetric waiting time distributions; 3) strong 
positive correlations between waiting times of successive 
events; and 4) weak or no correlation between intensities 
and waiting times. These statistical similarities, together 
with the fact that the crustal energy liberated by star- 
quakes is sufficient in principle to fuel the soft gamma- 
ray flashes, suggest that SGRs are indeed powered by 
starquakes. Moreover, there is also a strong correlation 
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between the magnitude and waiting times both for ac- 
tive earthquake regions and for SGR 1806-20. Thus, 
the statistical similarities between earthquakes and SGR 
events argue for physically similar origins. The Quasi- 
Periodic Oscillations (QPOs) observed at late-times of 
giant flar es of SGR 1806-20, SGR 1900+14, and SGR 
0525-66 (jMareghettil [200l iWatts fc Strohmaverl 12001 
constitute another piece of observational evidence in fa- 
vor of starquakes. These QPOs are most likely due to 
seismic oscillations induced by the large crustal fractures 
occurring in extremely energetic events similar to what 
happens after earthquakes. Such oscillations could be 
limited to the cru st or involve the ent ire neutron star. 

Crackling noise iSethna et a l. (2001) arises when a sys- 
tem responds to changing external conditions through 
discrete, impulse events spanning a broad range of sizes. 
Bak, Tang, and Wiesenfeld (1988) introduced a connec- 
tion between dynamical critical phenomena and crack- 
ling noise. They emphasized how systems may end up 
naturally at the critical point through a process of self- 
organized criticality Based upon this idea of the crack- 
ling noise, Kondratyev (2002) studied the statistics of 
magnetic noise in neutron star crusts, and compared its 
intensity and statistical properties to the burst activity 
of SGRs. He then argued that the noise could originate 
from magnetic avalanches. However, because of the re- 
quired inhomogeneous crust structure, he postulated the 
existence of magnetic domains within the neutron-star 
crust for an interior magnetic field strength in the range 
~ 10 16 — 10 17 G. Using the randomly jumping interact- 
ing moment (RJIM) model, it was shown that the burst 
intensity and waiting time distributions are not only in 
good agreement with observations, but also are analo- 
gous with the statistical properties of SGRs. 

Whether or not magnetars are the source of SGRs and 
AXPs, as relics of stellar interiors, the study of the mag- 
netic fields in and around degenerate stars should give 
important information on the role such fields play in star 



Suh & Mathews 



formation and stellar evolution (Suh & Mathews 2001a). 
The scalar virial theorem implies an allowed internal field 
strength for a star of B < 2 x 10 8 (M/M Q )(R/R Q )- 2 G 
for a star of size R and mass M. For a typical neutron 
star the maximum interior field strength could thus reach 
B ~ 10 18 G. Since strong interior magnetic fields modify 
the nuclear equation of state f or degenerate stars, their 
structure will also be changed ()Cardall et al.fe OOl). 

Even though SGRs appear to be observable conse- 
quences of starquakes or surface fractures, their detailed 
mechanism is not known. Moreover, it is unlikely that 
AXPs could only be explained by starquakes in strong 
magnetic fields. Therefore, in this paper we suggest a 
magnetic domain model to correlate smoothly between 
the statistics of starquakes and magnetic avalanches in 
magnetar crusts. In this work magnetic properties of 
magnetar-matter such as the magnetization and the sus- 
ceptibility are calculated in the framework of three dif- 
ferent representative equations of state. We consider 
an ideal n — p — e gas, relativistic Hartree mean field 
theory, and the ma gnetic Baym, P ethick, and Suther- 
land ( BPS) model (lLai fe Shapiro! Il991h. It has been 
shown (jBroderick et al.l 2000) that the magnetization of 
magnetar-matter undergoes large de Haas — van Alphen 
oscillations. The magnetic susceptibility can then lead to 
a region unstable to the formation of magnetic domains. 
It has not yet been demonstrated, however, that mag- 
netic domains actually form in magnetar-matter. Here 
we show that it is indeed possible to form such magnetic 
domains in magnetars, and that these could affect the 
surface properties and structure of magnetars possibly 
leading to observable consequences such as starquakes, 
glitches, and X- or 7-ray emission. 

2. THE DIFFERENTIAL SUSCEPTIBILITIES FOR A 
MAGNETIC NPE GAS 

The magnetic equation of state for magnetar-matter 
has been described in Suh & Mathews (2001a) for an 
ideal npe gas. The magnetization of simple magnetar- 
matter m aterial can be derived from th e thermodynamic 
potential (Blandford & Hcrnauistl ll982l ). Brodricket. al. 
(2000) have generalized the formalism for a multicompo- 
nent system including interacting nucleons. Hereafter 
we introduce the following notation for magnetic fields. 
For material in a uniform magnetic field, the magnetic 
field H is relat ed to the flux density B by the relation 
(TPinpardlUQSOh : 



B = H + 4tt{1 - V)M(B) 



(1) 



where T> is the demagnetization coefficient that is fixed 
by the geometry of the system. For example, P«0 for a 
neutron star crust permeated by an approximately ver- 
tical magnetic field. In chemical equilibrium, the total 
magnetization M. is given by a simple sum of the con- 
stituent magnetizations. 



M= J2 M > 



(2) 



The magnetic susceptibilities are then given by \ = 
M./H, and the differential s usceptibility n is defined 
(jBlandford fc HernquistllrMl by 
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where fi is the chemical potential, T is the temperature, 
and V is the volume of the system. The total differen- 
tial susceptibility for magnetar-matter above the neutron 
drip density is then given by 



i] = 2^ Vj> where j = n,p, e 



(4) 



For a cold ideal npe gas, we can obtain simple expres- 
sions for the differential susceptibilities of the various 
components. For the electron differential susceptibility 
rje we have, 
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where fh e = ,/l + 2^ e n e r. For protons, 
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where A p = nK / J 1 + 2"/ p n P f — s p z k p and 

Jl + 2 lp n p f . In Eqs. (5)- (6), n f ' p = r 

where n = 1,2,3,... denotes the Landau levels and 
s e z ,p are the electron (e) and proton (p) spin projection 
on the magnetic field direction, j e , p — B/Bp p where 
Bp p = eh/m 2 c 3 are the quantum critical field for elec- 
trons and protons, and e ep — Ep P /m e p c 2 with the elec- 
tron and proton Fermi energy Ep P , respectively. 
Finally, for neutrons we obtain, 
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where 
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with m n — H-2s"k„7„ and 7„ = B/B™, B™ = eK/1 
In Eqs. (5)-(7), a — e 2 /he is the fine structure constant, 
s™ is the neutron spin projection in the magnetic field di- 
rection, and k p , K n are the anomalous magnetic moments 
for protons and neutrons respectively, as given below in 
Eq. (10). Here we use the same notation as in Suh & 
Mathews (2001a) for the particle Fermi energy and mag- 
netic field strength. 
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3. THE DIFFERENTIAL SUSCEPTIBILITIES IN THE 
RELATIVISTIC HARTREE THEORY 

For a system of strongly interacting baryons (neutrons 
and protons), the relativistic mean field (Hartree) theory 
should be a reasonable approximation for the descrip- 
tion of t he equation of state fo r magnetar-matter at high 
density (jBroderick et all 120001 IChakrabartv et alj|1997l) 
through the exchange of a and vector u, p mesons in a 
strong magnetic field. In the baryon Lagrangian for the 
relativistic Hartree theory, the anomalous magnetic mo- 
ments are included through the coupling of the baryons 
to the electromagnetic field tensor with a^ — | [7^1 7i/] 
and the strengths n p and n n given by 



9 P 



2m p c I 2 



ffn 



2m n c 2 



(10) 



where g p = 5.58 and g n — —3.82 are the Lande ^-factors 
for protons and neutrons, respectively. In this work, we 
can ignore the possible scalar a, the vector u) and the iso- 
vector p meson self-interactions. Therefore, although the 
electromagnetic field is included in the total Lagrangian, 
it assumed to be externally generated (and thus has no 
associated field equation) and only frozen-field configu- 
rations will be considered. The effective baryon mass 
f^b=p,n is then given by the coupling to the a meson, 



it = m b -( 9a /m a f(n s p +n s n ) , 



(11) 



where g a and m a are the a meson coupling constant and 
mass respectively. In Eq. (11), n p is the scalar number 
density for protons, 
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(m e /m*) 2 7e, while the scalar number density for neu- 
trons is 
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with m* = 1 + s"K n 7n* and 7* = (m e /m*) 2 7 e . For sim- 
plicity, the nucleon rest mass is taken as m = m n = m D 
in the numerical calculation (Chakrabar tv et al.lll997ft . 

Assuming a mixture of neutrons, protons, and elec- 
trons in chemical equilibrium, the chemical potentials 
are related by 

Pn = Pp + Pe , (14) 



while the condition of charge neutrality gives 
n p (e p ,j p ) = n e (e e ,7e) ■ 



(15) 
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Fig. 1. — Variation of the effective baryon mass m* /m in the 
relativistic Hartree model as a function of baryon density for 
magnetic field strengths of f e = 0.01 (solid line) and 10 (dash 
line) . 



Given the nucleon-meson coupling constant and the co- 
efficients in the scalar self-interactions, the field equa- 
tions can be solved self-consistently for the chemical po- 
tentials, pj (j = n,p,e), and the meson field strengths 
in a uniform magnetic field B along the z axis corre- 
spondin g to the choice of the g auge for the vector poten- 
tial A^ (jBroderick et a.1.1 l2000f ). In this work, we adopt 
the following coupling constants and mesons masses: 
gl{m N /maf = 357 47, gl(m N /mu>) 2 = 2 73.78, and 
g 2 p (m N /mp) 2 = 97.0 (|fiorowitz k. Serotj|1981[) . 

Figure 1 shows the effective baryon mass m* /m as a 
function of baryon density p for magnetic field strengths 
of 7 e = 0.01 (solid line) and 10 5 (dash line) calculated in 
the model of Horowitz and Serot (1981). For a magnetic 
field strength less than ~ 10 18 G, This figure shows that 
the effective nucleon mass is not significantly affected by 
magnetic field strength. Broderick et al. (2000) and 
Chakrabarty et al. (1997) have obtained similar results. 
This effective baryon mass modifies the baryon dispersion 
relation in dense magnetar-matter. 

4. THE DIFFERENTIAL SUSCEPTIBILITY IN THE 
MAGNETIC BPS MODEL 

For matter in thermodynamical equilibrium below the 
neutron drip density, pdrip ~ 4.3 x 10 11 g/cm 3 , we adopt 
the magn etic Baym , Pethick, and Sutherland (BPS) 
model (lLai fc Shapiro! 1 1991D and use the semi-empirical 
mass formula (Shapiro & Teukolsky 1983). For simplic- 
ity, we only consider |gFe nuclei in the numerical cal- 
culation. Then, the magnetization and the differential 
susceptibility for the magnetic BPS equation of state are 



Suh k Mathews 




lO" 1 10° 


10' 


io 2 

B/B' 


io 3 


io 4 


10"" 



10 



Fig. 2. — Region of the magnetic domain formation instability 
for an ideal npe gas as a function of baryon density for various 
magnetic fields in units of B^ = 4.414 X 10 14 G. The region above 
the solid line denotes the parameter values for which conditions 
are unstable for the formation of magnetic domains. The dash line 
indicates the nuclear saturation density pnuc ~ 2.8 X 10 14 g/cm 3 . 



given by 

Mbps = M e 



-Mi 



rjBPS = Ve 



;Vl 



(16) 



where M. e is the magnetization of the electron gas and r\ e 
is given in Eq. (5). In Eq. (16), A4l is the magnetization 
for the bec Coulomb lattice energy. Then we drive here 
the lattice differential susceptibility tjl to be, 
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(17) 
where Z is the average atomic number of the nuclei. 

5. MAGNETIC DOMAIN FORMATION 

In general, the magnetization of a system is small com- 
pared with the external magnetic field H. However, when 
the system is sufficiently cool so that its thermal energy 
is smaller than the spacing of the Landau levels, the mag- 
netization can undergo large de Haas— van Alphen oscil- 
lations with either changing magnetic fields or a changing 
Fermi energy. Under these conditions, it sometimes be- 
comes energetically favorable for the system to separate 
into two phases containing differe nt flux densitie s. This 
is the so-called Schoenberg effect (|Pippardl[l9 80). This 
means that although -lirx is less than unity for a mag- 
netized npe gas, in certain regions 47777 can exceed unity 



which implies the possible exis tence of magnetic domains 
(|Blandford k Hernquistlll982D . That is, when the differ- 
ential susceptibility obeys n > 1/47T, and dH/dB < 0. 
Then magnetar matter in thermodynamic equilibrium 
becomes unstable to the formation of magnetic domains 
of alternating magnetization. For the case of a vanish- 
ing demagnification coefficient, T> = 0, the material will 
separate into two phases corresponding to different mag- 
netization. 

For magnetar-matter above the neutron drip density 
Pdrip, we considered an ideal pure non-interacting cold 
npe gas as well as the relativistic Hartree model. How- 
ever, in the density region between pdrip and p nuc ~ 
2.8 x 10 14 g/cm 3 , neutron-star matter is composed elec- 
trons, nuclei, and free neutron gas so that we can not 
directly apply the ideal npe gas model in this density 
regime. For example, for non-magnetic neutron-star 
matter in this intermediate density regime, we can em- 
ploy t he Baym, Bethe, and Pethick (BBP) equation of 
state ()Bavm et al.| [l988') . This BBP model is based upon 
a compressed liquid drip model of nuclei. It gives some 
corrections to the ideal npe equation of state. There- 
fore, if we adopt the magnetic BBP model, the region 
below the dash line in Figs. 2 and 3 will be shifted 
to the left because of the BBP equation of state [See 
I Shapiro k Teukolskvl ([1983)] . This means that for a fixed 
magnetic field strength the density region in which mag- 
netic domains can be formed increases. However, since 
there is no physical model in the intermediate density 
regime with a strong magnetic field, we can simply de- 
scribe this regime using an analogy from our ideal npe 
gas model. 

Figure 2 shows the regions of matter density and mag- 
netic flux density for an ideal npe gas above the pdrip 
in which magnetar-matter is unstable to a phase separa- 
tion into magnetic domains. For magnetic field strength 
of ~ 10 14 — 10 15 G, magnetic domains cannot be formed 
in the lower density region below ~ IO 13 g/cm 3 . How- 
ever, it is possible for magnetic domains to form in the 
density region higher than ~ IO 13 g/cm 3 . 

Figure 3 shows the unstable region for a relativistic 
Hartree mean field model in which the baryon effective 
mass is taken into account. When we consider the effec- 
tive baryon mass within the relativistic Hartree theory, 
domain formation can significantly occur above a den- 
sity of p ~ 10 14 g/cm 3 and a field strength B ~ 10 16 G. 
The effective baryon mass lowers the density at which 
magnetic domain formation occurs in the core of a mag- 
netar in which strong magnetic fields are expected. We 
also find that magnetic domain formation could not be 
formed for a magnetic field strength of < IO 18 G. 

Finally, in Figure 4, the magnetic domain instability 
regions obtained by using the magnetic BPS model are 
depicted in the outer crust of magnetars. Unlike the 
other two cases, in this case the magnetization is domi- 
nated by electrons. Therefore the conditions for the on- 
set of the magnetic domain instability is affected by the 
occupation of Landau levels. This causes the instability 
conditions to vary as a function of field strength for fixed 
density. In Figure 4 we see that magnetic domains can- 
not be formed if the density is less than p ~ 5 x 10 10 
g/cm 3 for a magnetar having a typical surface magnetic 
field of B ~ 10 14 ~ 15 G. However, in the density region 
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Fig. 3. — Region of the magnetic domain formation instability 
for the relativistic Hartree theory in which the effective baryon 
mass is taken into account. The region above the solid line 
denotes the parameters for which the conditions are unstable for 
the formation of magnetic domains. The dash line indicates the 
nuclear saturation density p n uc ~ 2.8 X 10 14 g/cm 3 . 
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Fig. 4. — Magnetic domain formation instability for IgFe in the 
magnetic BPS model. The shaded regions denote the parameter 
regions for which conditions are unstable for the formation of 
magnetic domains. 



around the neutron drip density, magnetic domains can 
be formed. We can also see magnetic domain formation 
in regions at low density > 10 11 g/cm 3 and at a rela- 
tively low magnetic field strength > 10 13 G. This is the 
case considered by Blandford and Hernquist (1982). For 
magnetic white dwarfs, Adam (1986) also obtained the 
similar result for a magnetized electron gas. 

We thus find that possible unstable regions for mag- 
netic domain formation are within the deep outer crust, 
deeper part of the inner crust, and core of magnetars. 
This means that the outer shell at low density less than 
p ~ 10 10 g/cm 3 is stable against magnetic domain forma- 
tion and should consists of strongly magnetized material 
without magnetic domains. 

6. SGR AND AXP MECHANISMS 

This magnetic domain formation might be an impor- 
tant clue to explain SGRs and AXPs in the magnetar 
model. As the density increases, the number of Landau 
orbital also increases. Therefore, Based upon Blandford 
& Hermquist (1982), we can calculate the thickness of 
the crust, Az, using the gravitational potential energy 
of a nucleon and the Fermi energy e e per electron. We 
also know the maximum number of the Landau orbitals 
n = (el — l)/27 e , where 7 e = BjB e c . For a fixed B, 
An ~ e e Ae e /B. Therefore, we can get Az/An. For 
An = 1, we finally can estimate the spacing between 
layers associated with the maximum Landau level: 




(18) 
where g is the surface gravity, /x e is the mean molecu- 
lar weight per electron, and e e = Ep/m e c 2 denotes the 
electron Fermi energy. Under the conditions for domain 
formation, the spacing Az would be the average vertical 
scale between domain interfaces as well as a horizontal 
scale of the domains if they are in local pressure equilib- 
rium. [However, the actual size and shape of the domains 
are difficult to determine.] Now, for a fixed the Fermi en- 
ergy, we get an equation An which is given by AB. For 
An = 1 then the horizontal variation in the magnetic 
field will be 
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(19) 



This implies that there are different magnetic flux den- 
sities in region where magnetic domains can form. How- 
ever, in the outer crust at lower density, the magnetic 
fields are homogeneously distributed and tightly pinned 
to the matter. 

Ac cording to the magnetic domain theory (jPippardl 
Il980h . the magnetic domain walls move or grow, and the 
magnetic domains rotate within the material. There- 
fore, the movement and rotation of magnetic domains 
can cause a physical dimensional change and produces 
the maximum possible strain on the crust material. This 
process is very similar to that of the internal structure 
of the earth in which a sudden collapse or strain of the 
mantle below the earth's crust sometimes occurs. 

In regions where magnetic field distortion increases, 
magnetic domains could be formed. At the boundary 
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there will also exist regions around each wall where the 
magnetic field is distorted. The formation or adjustment 
of domain structure involves magnetic field fluctuations 
of a few percent amplitude which have an anisotropic 
magnetostrictive stress 2ttA4 2 associated with the mag- 
netization. Any sudden readjustment of the domain 
structure will cause a local departure from isostasy 
which will be relieved on an oh mic dissipation timescale 
(rn ~ (tA/Attc 2 ~ 10 4 yr) (jBlandford fc Hernquistl 
1982). These anisotropic magnetostrictive s tresses may 
be la rge enough to crack the outer crust (Blaes et al.l 
QUI). 

Then, we can estimate the physical length variation 
of the magnetic domains. The bulk modulus (Ks) is 
defined as the pressure increase needed to effect a given 
relative decrease in volume. For a gas, the adiabatic bulk 
modulus Ks is approximately given by K$ = TP, where 
r is the adiabatic index. Then, 
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for r ~ 1.8 and x — 1 at the neutron drip density (4.4 x 
10 11 g/cm 3 in the magnetic BPS model (Lai & Shapiro 
119911 )). Therefore, with A « 0.6 for Fc (Stc wartl 11951 )7 
the length change due to magnetic domain formation is 
finally given by 



^ = 1^,0.02 
I 3 V 



(21) 



This means there is a 2 m length change for a 100 m 
characteristic domain size when a magnetic domain is 
formed in the deep outer crust of a magnetar. Now we 
can estimate a cracking timescale in the outer crust to 
be tc ~ Az/v s ~ 0.1 ms, where v s — y /Y/p (with 
Y = shear modulus) is the shear velocity (Blacs et al. 
1989). Kontratyev (2002) has analogized this cracking 
timescale as the avalanche spanning time which is con- 
sistent with the rise time for SGR giant burst s (jHurlevi 
2000J iKouveliotiu et~aT1[200l lMareghettill200§ ) . 

Finally, we can estimate the elastic energy released, 
AEjj ~ Az R\)<j, where Ro ~ Az is the characteristic 
horizontal size of the domain, using the magnetic stress 
energy a ~ X 2 B 2 - Then, we obtain the released elastic 
energy 



AE D - 6 x 10 4 V<9„ 
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10 4 cn 
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erg. (22) 



the crust by magnetostrictive stress would be the mech- 
anism of the observed SGRs. We suggest that magnetic 
domain formation and any sudden readjustment of the 
domains can produce an energy source for soft gamma- 
rays in SGRs and X-rays in AXPs. 

However, there is evidence that AXPs have stronger 
magnetic fields than SGRs (jMareghettil I2008D . With 
stronger magnetic fields, it would be hard for the mag- 
netic domains to be formed in the outer crust of mag- 
netar. This means that there would be little cracking 
of the outer crust by the magnetostrictive stress. The 
possibility remains, however, that AXPs could produce 
a giant blast like SGRs in this magnetic domain model. 

7. SUMMARY 

In this work we have studied magnetic domain forma- 
tion as a new mechanism for SGRs and AXPs. In this 
paradigm magnetar-matter separates into two phases 
containing different flux densities. We have identified 
the parameter space in matter density and magnetic field 
strength at which there is an instability for magnetic 
domain formation and have shown that such instabili- 
ties are likely to occur in the deep outer crust for the 
magnetic BPS model, and in the deeper part of the in- 
ner crust and core for magnetars described by relativis- 
tic Hartree theory. Moreover, we have estimated the 
strain on the outer crust induced by the formation of 
such domains and found that the anticipated energy re- 
lease is comparable with the energy emitted by typical 
SGRs. Hence, we propose that the magnetic domain 
formation scenario described here represents a new pos- 
sible mechanism to drive the giant flares of SGRs as well 
as X-ray outbursts and the quiescent phase of AXPs. 
At the very least, this proposal warrants further inves- 
tigation. Moreover, since the physical length variation 
caused by the magnetic domain formation might lead 
to solid crustal deformation and catastropic cracking, 
SGRs might be sou rces of gravitational waves (GWs) 
(|Abbott et al.l |2008[ ) even though there is not yet ev- 
idence of GWs associated with observed SGR bursts. 
However, if it becomes possible to detect GWs from 
SGRs, that may be a way to verify the magnetic domain 
model in magnetars. Clearly, the next step is to under- 
take detailed dynamical numerical studies of the forma- 
tion and evolution of such magnetic domains in neutron 
star crusts. Effor ts along this line are currently underway 
(jSuh et al.H20"09l) . 



where 9 m the maximum allowed strain angle, and \ = 
M./H is the magnetic susceptibility. This is also the 
typical energy released in SGRs. Hence, this cracking of 
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